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Fluctuations in the return time statistics of a dynamical system can be described by a new spec- 
trum of dimensions. Comparison with the usual multifractal analysis of measures is presented, and 
difference between the two corresponding sets of dimensions is established. Theoretical analysis and 
numerical examples of dynamical systems in the class of Iterated Functions are presented. 



PACS numbers: 05.45.-a 05.45.Df 05.45.Tp 

Suppose that T is a transformation of the space X into 
itself, which preserves the probability measure /x. Let A 
be a subset of X, and x a point in A. We let ta(x) be 
the (integer) time of the first return of x in A: 



t a (x) = mf{n > s.t. T n (x) £ .4}. 



(1) 



Poincare theorem guarantees that the return time is 
almost certainly finite with respect to any invariant mea- 
sure. It can be rather long, and particularly short as well: 
if x is a fixed point of T, then obviously ta(%) = 1 for 
any set A containing x. In this paper, we shall prove 
that the distribution of return times is characterized by 
multi-fractal properties which can be properly described 
by the tools of the so-called thermodynamical formalism 

Much recent research has been focused on the local 
statistics of these returns, in the following sense: Take 
for A a ball of variable radius e centered at the point x: 
A = B £ (x). Next, consider the cumulative distribution 
of the first return times, in the ball, of all points of this 
latter: Let A(e) = fj,(B e (x)), and define m(x,e,s) := 
H{y e B e (x) s.t. \(e)T Bc ( x )(y) > S }/-Me)- In man Y in- 
stances, these statistics have an exponential character: 
m[x,e, s) — > e~ s , as e tends to zero. Rigorous proofs of 
this fact have been produced under an ever lessening set 
of hypotheses || . 

More fundamental for our analysis than the previous, 
is Kac theorem a classical result of the local analysis: 
it predicts that, whenever the measure [i is ergodic with 
respect to T, for any measurable set A (and not just 
balls), the expectation of ta over the set A is just the 
inverse of the measure of A: 



TA{y)—r-rr = — r-TT- 

A n{A) ll{A) 



Furthermore, one must quote the theorem of Ornstein 
and Weiss B: Let A be a generating partition of X, and 
let us refine it around all the points x G X by defin- 
ing the cylinder of order n, A n (x), as the intersection 
of all the elements of A, T^ 1 A,. . . ,T~ n+1 A containing x. 
The Ornstein- Weiss theorem states that, whenever the 
measure jj, is ergodic, the following limit exists /i-almost 
everywhere, and is equal to h([x), the metric entropy of 



lim 



(3) 



This remarkable result is historically the first link be- 
tween a thermodynamic quantity, and the return times. 
Parallel to the Ornstein- Weiss theorem, the Shannon- 
Mc Millan-Breiman theorem gives the metric entropy by 
means of the exponential decay of the measure of cylin- 
ders around almost all points. It is common in the phys- 
ical usage to replace cylinders with balls: this allowed in 
force of the Brin-Katok theorem ||. The advantage of 
balls over cylinders, when considering practical applica- 
tions, is apparent. 

Let us therefore replace A n (x) and n in eq. (||) by 
B e (x) and — loge, respectively. In the case of Gibbs mea- 
sures of Axiom- A diffeomorphisms || , and of a wide class 
of maps of the interval [[7| one has been able to prove rig- 
orously that the modified limit (^) exists /z-almost every- 
where, and is equal to the information dimension _D M (1), 
of ji: 



lim 



\og(T Beix) (x)) 
loge 



(4) 



(2) 



The definition of _D M (1), and an informal proof of eq. (Q), 
will be given momentarily. 
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So far, the analysis has been local. Quite different-and 
more complex-is the case of the global statistics that we 
consider in this paper. We let balls of fixed radius e cover 
X, and we compute the integrals 

T T (e,q):= T B ~^ x) (x)dfi(x), (5) 

where q is a real quantity. These are the statistical mo- 
ments of the time required to come back to a neighbor- 
hood of the starting point: it is clear that local analysis 
alone has little to say about the scaling, in the radius e, 
of these quantities. 

Indeed, L T (e, q) are a sort of partition function, quite 
akin to those employed in the thermodynamical formal- 
ism. Their scaling for small e defines a set of dimensions, 
D T (q): 

T T (e,q)~e D ^~ 1 \ (6) 

that we call return time dimensions. The study of these 
dimensions is the object of this Letter. In level of impor- 
tance, the first question is about their relations with the 
usual quantities of the thermodynamical formalism. 

Let us first observe that the return time T Be ^(x) can 
be interpreted as a "single-point" sample of the integral 
in Kac theorem, eq. (^). Therefore, one could estimate 
that 

T Be{x) {x)~ »{B e {x))-\ (7) 

It is worth recalling now that the local dimension a(x) 
of the measure /i at the point x is defined by the scaling 
fj,(B E (x)) ~ e a ( x \ in the limit of small e. Moreover, the 
relation a(x) = D^{\) holds almost everywhere in X, as 
discussed above. This fact, and eq. (0), then imply the 
theorem expressed by eq. (||). 

The approximate equality (^) has been already 
adopted in to evaluate via eq. (||), the exact ther- 
modynamical partition function T^ (£,<?), 

T^s,q):= f / /(B £ (x))"- 1 d/i(a:)~ £ ^W^- 1 ), (8) 
Jx 

whose scaling for small e gives the well-known spectrum 
of generalized dimensions D^{q) |l(| . 

The substitution of eqs. (||) by eqs. ^,0) is partic- 
ularly significant, for it allows one to treat dynamical 
systems endowed with physical measures of the Sinai- 
Bowen-Ruelle type: in this case the integral can be re- 
placed by a Birkhoff sum over the trajectory x\ := T 1 (xq), 
I = 0, . . ., of a generic point xq\ 

n — l 

T T (e,q)= lim t£ ( * (9) 

Indeed, this Birkhoff procedure, eq. @, is the one orig- 
inally employed in where it is claimed to produce 
the spectrum of measure dimensions D^(q). 



It must be remarked that eq. (|9|) considers a single sum 
along the orbit, while in the usual Grassberger-Procaccia 
technique for the g-correlation integrals, one must work 
with a double summation, where all the couples of points 
along the same (or different) orbit are compared O]. 

We are now ready to introduce the main result of this 
paper: contrary to the usag e of |,§, the estimate (@) 
does hold only for the first moment, as in Kac theorem, 
but not in the stronger sense required in eqs. (|B|), and 
(||): the variables Tg e / x -\ and ^,(B £ (x))~ 1 have different 
large deviation properties. Therefore, D^{q) and D T (q) 
are not the same function, and the latter defines a new 
bona fide spectrum of dimensions. 

The abstract proof of this statement would be little 
informative, if not paralleled by a specific example. It is 
therefore convenient to introduce a family of dynamical 
systems whose invariant measures are completely known, 
in the sense that the spectrum of dimensions D^q) can 
be easily and precisely computed: these are the so-called 
systems of iterated functions, or I.F.S., see [Q for de- 
tails. In the simplest, one-dimensional case, an I.F.S. 
consists of a collection of contracting maps of the line, 
(f>i : R — > R, i = 1, . . . , M, where M is the number 
of maps. These latter can be thought of as the inverse 
branches of the dynamics: (T o <j)i){x) = x for any i. In 
so doing, T is characterized by a mixing repeller, that is 
also the attractor of the collective action (in the sense 
that will be made soon clear) of the maps in the I.F.S.. 

A family of invariant measures for T can be con- 
structed assigning a probability, iXi, to each map <pi, 
^2iLi n i = 1- The dynamics of I.F.S. can be constructed 
by sequentially applying maps <p a , where a is chosen in 
a random fashion with probability ir a . Any orbit of the 
probabilistic I.F.S., when time-reversed, becomes an or- 
bit of the deterministic dynamics T . Notice that return 
times are invariant under time reversal. 

The reader will find useful to consider the usual ternary 
Cantor set measure as the invariant measure of the fol- 
lowing 2-maps I.F.S.: 4>i(x) — x/3, 4>2(x) = (x + 2)/3, 
tti = 7T2 = 1/2. The related invariant measure is an 
example of mono-fractal: all generalized dimensions are 
equal to the constant S = j§|f ■ Let us perform the re- 
turn time analysis of this dynamical system following the 
Birkhoff sums approach, eq. (^). Other techniques of es- 
timating the integral (||) will be presented elsewhere. 

We can first verify that the scaling relation (||) holds: 
fig. ^] assures us that this is indeed the case. 

We can then extract the return dimensions D T (q) ver- 
sus q: the numerical results are plotted in fig. ||. They 
reveal a non-trivial range of dimensions: as far as return 
times are concerned, the dynamics is multifractal. More- 
over, the return time dimensions are consistent with the 
constant value S (within the statistical error bars) for 
negative dimensions, but are significantly different from 
this latter for positive values of q, the more, the larger 
the value of q. 

Resorting now to exact analysis, we can confirm these 
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results: precisely, we are able to prove that i) -D M (0) = 
D T (Q), and ii) D T (q) — > for large q. The second fact 
proves that the return dimensions D T (q) are significantly 
different from measure dimensions, while the first hints 
at relations that, at least in certain dynamical systems, 
may exist between them. 

To see that i) is true, observe that the repeller of T 
can be hierarchically covered by intervals Af, of length 
l/3 fc , where k = 0, . . . is the index of the generation, and 
I = 1, . . . , 2 k is the label of the single interval within its 
generation. This is simply the usual construction of the 
ternary Cantor set, in which one carves a hole of length 
l/3 k in the middle of each interval generation k — 1. 

For a fixed value of k, let us choose e = £k = l/3 fe . 
Then, simple geometric considerations assure us that 
whenever x belongs to A k , with k and I fixed, it is also 
contained in the ball of radius Sk centered around any 
other point in Af. In addition, these balls intersect no 
interval of generation k other than Af. This permits us 
to employ Kac theorem, eq. (S), to compute exactly the 
sum r T (e fe ,0): 



r r (e fc ,o) 



1" 

E 

i=o 



T A h{x)dn{x) 



(10) 



It is then clear that i) holds. 

Remark that the intervals Af , are also the cylinders of 
order k of the dynamics. Therefore, the result we have 
just proven also applies to the return time dimensions, 
when defined via cylinders, and predicts that the dimen- 
sion of order zero is equal to the topological entropy. 
Following the same idea, one also expects that the cylin- 
der return time dimensions will be related to the Rcnyi 
entropies of the measure [|l3| . 

Observe that the above results are by no means re- 
stricted to the one-dimensional case. At the same time 
it is evident that use of Kac theorem is permitted only 
for q — 0, for otherwise different fluctuation properties of 
the return time statistics set in. This leads us directly to 
the statement ii). 

Remark that the return times define partitions of X by 
the sets Rn{s) := {x G X s.t. tb s ( x ){x) = n}. This allows 
us to re- write the moments r r (e, q) as a new summation: 



r r (e,?) 



n=l 



n L -ifi(R n (e)). 



(11) 



When q is large, terms with small n lead the sum (11). 
Let us retain just the first of these: it does not vanish 
by increasing q. This term is the measure of R\(e), the 
set of points that move less than e in a single iteration. 
Clearly, all fixed points of T belong to this set. Let x 
be any one of these, and let us assume that e is small 
enough, and that the transformation T is smooth. Then, 
Ri(e) contains the ball of radius e/A centered at 5, where 
A is the largest singular value of the matrix 1 — T', 1 is 
the identity, and T' is the Jacobian matrix of T at x. 



Moreover, the ball of radius e/a, where a is now the 
smallest singular value, contains the connected part of 
Ri(e) around x. Then, we can estimate that, for large q, 



T T (e,q)~ »{ B pe{x)), 



(12) 



x S.t. T(x)=i 



where the sum is extended to all fixed points of T, and 
where p is a multiplicative factor, which depends on x: 
a < p < A. When the number of terms in the sum 
is finite, the leading contribution to (|l2|) comes from the 
fixed point x with the smallest local dimension, a m : since 
p,(B e (x)) ~ e a ™, it follows at once that 



D T (q) 



9-1 



when q — > oo. 



(13) 



We can readily verify that eq. ( |13[) is satisfied by the 
Cantor set dynamics introduced above. The fixed points 
of T are here zero, and one, and they are both character- 
ized by the same value of the local dimension, a m = . 
In fig. U the relation (|l|) is validated by the numerical 
results of D T (q) for large q. 

This investigation is not limited to a single dimen- 
sion. We now add a two-dimensional example with a 
non-trivial dimension function D^{q): the motion on 
a Sierpinsky gasket, corresponding to the I.F.S. maps 
M*,V) = (x/2,y/2), <h(x,y) = ((x + l)/2,(y + l)/2), 
and 4>z(x,y) = (x/2, (y + l)/2), with non-uniform prob- 
abilities 7Ti = 4/10, 7T2 = 7T3 = 3/10. In fig. H we have 
compared the exact thermodynamical function D^(q) 
and the numerically evaluated D T (q) . For negative values 
of q the two dimensions are very close, and discrepancies 
may be due to the finite length of the Birkhoff sum, and, 
more importantly, of the fitting interval in e, which can- 
not be accounted for by the statistical definition of the 
error bars. For positive values of q, to the contrary, the 
two dimensions are rather different, and the asymptotic 
formula ( p"3] ) soon becomes a very close approximation of 
D T (q). 

In conclusion, we have shown that the statistics of re- 
turn times are characterized by multi-fractal properties, 
well defined by a new spectrum of dimensions D T (q). We 
have also found that, contrary to implicit previous state- 
ments in the literature, this spectrum does not coincide 
with the usual multi- fractal spectrum D M (q). It is now 
interesting to understand the meaning of the Legendre 
transforms associated to these new dimensions. We have 
proven that D T (0) = D^(0) under certain hypotheses 
on the dynamical system. Approximate equality seems 
to hold for negative q's as well This might pro- 

vide a means for computing negative measure dimen- 
sions, which is known to be a challenging numerical task. 
Finally, we have shown that in a large class of systems 
the relation D T {q) ~ holds for large q, where a m is 
the smallest local dimension at the fixed points of T. 
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FIG. 1. Return time partition sums F T (e,q) versus e for 
a set of equally spaced values of q ranging from q = —19 
(bottom curve) to q = 19 (top). Lines in between points are 
solely intended to connect data with the same value of q. 



FIG. 2. Generalized dimensions D T (q) versus q for the 
ternary Cantor set dynamics. Error bars are defined via the 
uncertainty in the least square fits of Fig. (|l|), under the 
usual statistical assumptions. The horizontal line is the (flat) 
measure thermodynamics : -D M (g) = log 2/ log 3. The dashed 
curve for q > is eq. (jl3]) with a = log 2/ log 3. 
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FIG. 3. Generalized dimensions D tl (q) (dotted line) and 
D T (q) (error bars) versus q for the Sierpinsky dynamics de- 
scribed in the text. Data have been computed by least square 
fits of Birkhoff sums with n = 800000. Also reported for 
positive q is the asymptotic formula ([l3|). 



